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Abstract. WR 1 18 is a highly evolved Wolf-Rayet star of the WC10 subtype surrounded by a permanent dust shell absorbing 
and re-emitting in the infrared a considerable fraction of the stellar luminosity. We present the first diffraction-limited 2.13 /mi 
speckle interf erometric observations of WR 118 with 73 mas resolution. The speckle interferograms were obtained with the 6 m 
telescope at the Special Astrophysical Observatory. The two-dimensional visibility function of the object does not show any 
significant deviation from circular symmetry. The visibility curve declines towards the diffraction cut-off frequency to ~ 0.66 
and can be approximated by a linear function. Radiative transfer calculations have been carried out to model the spectral energy 
distribution, given in the range of 0.5-25 /im, and our 2.13 visibility function, assuming spherical symmetry of the dust 
shell. Both can be fitted with a model containing double-sized grains ("small" and "large") with the radii of a = 0.05 and 
0.38/im, and a mass fraction of the large grains greater than 65%. Alternatively, a good match can be obtained with the grain 
size distribution function n(a) ~ a -3 , with a ranging between 0.005/im and 0.6/im. At the inner boundary of the modelled 
dust shell (angular diameter Oi n = (17 ± 1) mas), the temperature of the smallest grains and the dust shell density are 1750 K 
± 100 K and (1 ± 0.2) • 10 -19 g/cm 3 , respectively. The dust formation rate is found to be (1.3 ± 0.5) • 10 -7 M0/yr, assuming 
Kvind = 1200 km/s. 
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1. Introduction 

The Wolf-Rayet star WR118 (= CRL2179 = IRAS 18289- 
1001) in the catalog of WR stars of van der Hucht (2001) be- 
longs to the latest subtypes of the carbon-rich sequence (WC). 
It is classified as WC9-10 (Allen et al. 1977, Massey & Conti 
1983). Its very prominent IR excess is attributed to the thermal 
emission of a warm carbon dust shell. No remarkable varia- 
tions in the dust emission have been registered over the last two 
decades (Williams et al. 1987, van der Hucht et al. 1996) and it 
has been designated as a "persistent" dust-maker (Williams & 
van der Hucht 2000). While the dust grains are momentum- 
coupled to the fast stellar wind, fresh dust grains should be 
steadily condensing at a definite distance from the star, i.e. the 
dust formation rate is generally constant on the time scale of 
several decades. 

Among the known galactic WR stars, WR 1 18 is one of the 
optically faintest objects, but one of the brightest ones in the 
infrared regime. This is naturally explained by its heavy inter- 
stellar reddening and the obscuration by its conspicuous dust 



shell. A strong interstellar 9.7 /im silicate feature is observed 
in its spectrum with an optical depth of rg.7 m w 0.7 (Roche 
& Aitken 1984, van der Hucht et al. 1996). Using the rela- 
tion Av/Tg.rum < 21.5 (Schutte et al. 1998), an extinction 
of Ay < 15 can be derived. Allen et al. (1977) and Williams et 
al. (1987) list Ay ~ 15 and ~ 12.8, respectively. In the optical 
region, the brightness of WR 1 18 is only roughly known. 

From spectrophotometry, Massey & Conti (1983) estimate 
the v narrow-band magnitude to be approximately 22, and 
Cohen & Vogel (1978) list a V broad-band magnitude of — 20. 
Williams et al. (1987) modelled the IR spectral energy distribu- 
tion (SED) of the WC stars, finding that they could be well fit- 
ted with an optically thin shell consisting of small (~ 0.01 /im) 
amorphous carbon grains with a density distribution of p oc 
r~ 2 . The absence of the silicon carbide feature at 11.3 /im in 
the spectra of WR stars exclude SiC to be a constituent of 
their dust shells. In this paper we present diffraction-limited 73 
mas speckle interferometric observations of the dust shell of 
WR 118. Radiative transfer calculations have been performed 
to model both the SED and the 2. 13 /im visibility. 
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2. Observations and data reduction 

The WR 118 speckle interferograms were obtained with the 
Russian 6 m telescope at the Special Astrophysical Observatory 
on September 27, 1999. The speckle data were recorded with 
HAWAII speckle camera (HgCdTe array, 256 2 pixels, sensitiv- 
ity from 1 to 2.5 /im, frame rate 2 frames/s) through an inter- 
ference filter with a central wavelength of 2. 13 /im and a band- 
width of 0.21 /im. Speckle interferograms of the unresolved 
star HIP 87540 were taken for the determination of the atmo- 
spheric speckle transfer function. The observational parame- 
ters were as follows: exposure time/frame 20 ms; number of 
frames 210 (130 of WR 118 and 80 of HIP 87540; 2.13 /im 
seeing (FWHM) l'/4; field of view 5'/l x 5'/l; pixel size 26.4 
mas. The resolution was 73 mas. The visibility function (mod- 
ulus of the Fourier transform of the object intensity) was de- 
termined with the speckle interferometry method (Labeyrie 
1970). Figure |] shows the two-dimensional 2. 13 /im visibility 
and the azimuthally averaged 2.13 fim visibility of WR 1 18 up 
to the diffraction cut-off frequency. There is no evidence for 
deviation from circular symmetry. The visibility curve can be 
approximated by a linear function and dropped to 0.66 at the 
diffraction cut-off frequency. 

3. The dust shell model 

Using the quasi-diffusion method of Leung (1975), as imple- 
mented in the CSDUST3 code (Egan et al. 1988), we have 
constructed spherically symmetric radiative transfer models to 
match both the observed SED and the 2.13 /im visibility. The 
input parameters of the model are: (i) the spectral shape and 
bolometric luminosity of the star (/(A)iboi), after passing the 
shell, /(A) iboi transforms to /'(A)Lboi; (ii) the inner shell ra- 
dius (r; n ); (iii) the relative thickness of the shell, i.e. the ratio 
of outer to inner shell radius (Y out = r out /ri n ); (iv) the den- 
sity distribution p(r); (v) the chemical composition and grain- 
size distribution; (vi) the total optical depth at a given reference 
wavelength. 

The input parameter of the radiative transfer model, which 
relates the stellar bolometric luminosity and the inner shell 
radius, is the incident flux at the inner boundary of the dust 
shell: F[ n — Lboi/(47rr? n ). It is based on the assumption that 
the temperature of the dust grains is controlled by the radia- 
tion field only. This means that if we change Lboi but keep 
Fi n = const (i.e. change r; n oc V-^boi, resp.) we obtain (i) 
the same SED; (ii) the same normalized surface brightness dis- 
tribution as a function of b/r- m (b/r m : impact parameter); and 
(iii) the same visibility curve as function of qQ[ n (q: spatial fre- 
quency; 8i n : angular diameter of the dust shell's inner bound- 
ary, b/2r- m = G/9j n ). For more details of self-similarity and 
scaling behaviour of infrared emission from radiatively heated 
dust, see also Ivezic & Elitzur (1997) and references therein. 
For the sake of clarity, we have incorporated Lboi and r m sepa- 
rately in the CSDUST3 code instead of using F ln only. During 
the SED fitting procedure, we kept Lboi constant., i.e. changed 
only r in . 

Keeping in mind that the distance is defined by the relation 

d = yj /'(A)L b oi/47rL obs (A) (L obs (A): observed (dereddened) 




J3 



> 

0.4 - 
0.2 - 

o.o I — 1 1 — ' ' 

2 4 6 8 10 12 

Cycles per arcsec 

Fig. 1. Two-dimensional 2.13 /im visibility (top) and az- 
imuthally averaged 2.13 /im visibility of WR118 (bottom) 
shown up to the diffraction limit. 



spectral energy flux), i.e. d oc y/Lboi, it can be shown that the 
angular surface brightness distribution of the model, in particu- 
lar 0; n = 2r- ln /d, does not depend on Lboi- In our calculations 
0; n is determined by comparing the modelled and observed K- 
band (A = 2.2 urn) flux. We note that, if the observed SED is 
properly fitted, we get the same value of Q ln conducting the 
flux normalization at any other wavelength. In principle, the 
bolometric flux can be used as well, but for WR 118 this ap- 
pears to be unfavorable due to the lack of optical and UV data. 

In our model calculations, we firstly assumed a bolomet- 
ric luminosity of Lboi = 4 • 1O 4 L . Luminosity-independent 
model parameters are given in Tab. [[} Removing the redden- 
ing from the observed J magnitude of WR 118 and taking the 
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Table 2. Predicted WR118 parameters: color V — J; dereddened and reddenend visual magnitudes Vq and V; distance d, 
bolometric luminosity Lboi; radius r m and density p m at inner dust shell rim; and dust mass-loss rate A/dust- 
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Table 1. Model dust-shell parameters. The models M 1, M2, 
Ml.l and M 2. 1 refer to double-sized grains with a = 0.05 pm 
and a = 0.38 /im, the model M 3 to a grain size distribution of 
n(a) oc a~ 3 with 0.005 /im< a < 0.6 /im. tv is the optical 
depth at 0.55 pm; r- m , T- m and 6i n are the radius, the tempera- 
ture (refering to the smallest grains) and the angular diameter 
of the inner dust shell rim, rao.05 and mo.os are the abundances 
of 0.05 pm grains by number and mass, resp., and Ak is the 
2.2 /im dust-shell albedo. 
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0.36 
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0.36 
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17 


0.40 


16.0 



model V — J color, we obtain the (unreddened) model Vq mag- 
nitude of WR 118. 

Adopting a WC9-10 absolute magnitude of My = —4.9 
(Williams et al. 1987), d and Lboi can be estimated. Taking this 
new value of Lboi (see Tab. ||) we recalculated the luminosity- 
dependent model parameters such as r; n , the rate of dust for- 
mation at the inner boundary of the dust shell (AfdustX and the 
density of the dust shell at the inner boundary (p; n ) using the 
relations r in oc \/Lbol> Mdust oc V^boi, and p m oc l/V^boi- 
Concerning the dust formation rate, an outflow velocity of 
V wind = 1200km/s was adopted (Williams et al. 1987). The 
model predictions for all these luminosity dependent parame- 
ters of WR 1 18 are given in Tab. |2|. 

For the SED of the star, /(A), we took basically a black 
body of 19000 K (van der Hucht et al. 1986). We corrected it 
in the region of 0.5-1.25 /im to have /(A) oc A -3 (Williams 
et al. 1998) corresponding to (V — J) « 0.0. In the infrared, 
corrections for the continuum radiation of the ionized gas shell 
have been made to obtain (J — K) w 0.25 and /(A) oc A~ 2 - 7 
longwards of 10 /im (Williams et al. 1987). It should be noted 
that these corrections have only minor effects on the model 
SED due to a very strong IR excess caused by the emission 
of dust grains. 

The density distribution of the dust shell is given by p oc 
r~ 2 , i.e by a uniform outflow model. The outer shell bound- 
ary was chosen to be Y out — 10 3 . Larger values of Y out af- 



fect the SED only in the far-infrared, i.e. far beyond the IRAS 
25 /im photometric band, which defines the long-wavelength 
tail in this study. Concerning the optical dust-grain constants, 
we considered the amorphous carbon grains (cellOOO) of Jager 
et al. (1998). The reference wavelength of the optical depth was 
chosen to be 0.55 pm, and the model SED is normalized to the 
observed one in the J band. 

Figure |2| shows the 2.13 pm model visibility for differ- 
ent optical depths (0.6 and 4) and grain sizes (single-sized 
grains with 0.01 pm, 0.1 pm and 0.38 pm) for a dust shell with 
Tin/R* = 200 in comparison with the observed visibility. Only 
for the model with the largest grains, the visibility curve ap- 
proaches to a linear function, i.e. its shape approaches the one 
of the observed visibility of WR 118. 




q [arcsec -1 ] 



Fig. 2. 2.13 pm model visibility function of WR 118 for dif- 
ferent optical depths and grain sizes for a dust-shell with 
r; n = 200i?* together the observed visibility. The dotted, 
short-dashed, dot-dashed, solid and long-dashed lines refer to 
t v = 0.6, a = 0.01/xm; t v = 0.6, a = 0.10/im; t v = 4.0, 
a = 0.10/im; Ty = 4.0, a = 0.38/im; and ry = 0.6, 
a = 0.38/im, resp. 

The transformation of the shape of the 2.13 pm visibil- 
ity curve from parabolic to rectilinear is connected with the 
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transformation of a mostly thermally emitting dust shell to a 
highly scattering one when the albedo of dust grains is get- 
ting large. The 2.13 /im albedo of the Jager et al. (1998) carbon 
grains with a=0.1 /mi is 0.04, i.e. the dust shell's radiation is 
dominated by direct thermal emission at this wavelength. The 
albedo of the corresponding 0.38 /im grains is 0.4, i.e. the dust 
shell scatters the IR light very efficiently. Increasing the opti- 
cal depth of an "emitting" and a "scattering" dust shell results 
into a slight increase and decrease of the visibility's curvature, 
respectively (Fig. ||). 

Keeping in mind the very hostile environment for dust con- 
densation in the shells of WR stars, it appears to be obvious 
to assume that the dust grains are rather small (Williams et al. 
1987). However, such models will not be able to reproduce the 
shape of the visibility curve of WR 118. On the other hand, it is 
clear that the dust shell cannot consist entirely of large grains, 
which grow up from their small progenitors. In our double- 
sized grains models, we considered the combination of small 
grains (a = 0.05 /im) and large grains (a = 0.38 /im). 

The temperature of the large grains is ~ 1.5 times lower 
than that of the small ones. Thus, if one increases the abun- 
dance of the large grains, r ln has to be decreased to preserve the 
SED fits in the near-infrared domain. Although this even wors- 
ens the condition of dust-grain condensation in the correspond- 
ing models, the observed shape of the visibility curve clearly 
requires to include large grains in the dust shell of WR 118. 

The decrease of the small grains' size does not affect the 
SED nor the visibility substantially, i.e. the grains with a < 
0.05 /tm behave similarly, in this respect, in the optical and IR 
domain. The increase of the relative abundance of the small 
grains mainly reflects in an increase of the emergent radia- 
tion's reddening in the optical region, and a decrease of the 
dust shell's albedo in the infrared. The latter leads to an es- 
sential change of the visibility's curvature in the "wrong" di- 
rection. Both these points define the abundance of small grains 
in the dust shell. However, in the case of WR 118, we have 
no accurate estimations of its brightness in the optical wave- 
length range, and thus cannot determine how much reddening 
we need. Calculations were performed adopting Ay = 12.8 
and 14, which leads to J « 4.65 and 4.31 (Allen et al. 1977, 
Williams et al. 1987), resp. 

The increase of the large grains' size does not affect the 
SED and the 2.13 /im visibility very much, either. However, 
the size of the large grains cannot be decreased significantly 
because the 2.13/im albedo of the Jager et al. (1998) grains 
with a = 0.25 /im is ~ 0.27 being already too low to allow a 
proper fit of the visibility. Besides double-sized grains models, 
we also calculated models with a power-law grain-size distri- 
bution of n(a) oc aT 1 (0.005 /tm< a < 0.6 /im) to investigate 
which exponent 7 is most appropriate. For all these models, 7] n 
is calculated for the smallest dust grains with a = 0.005 /tm. 
This temperature exceeded that of the a = 0.05 /im grains by 
- 5%. 

4. Results of calculations and discussion 

Figures show the observed 2.13 /zm visibility of WR 118 
and its SED corrected for interstellar extinction of Ay = 12.8 
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Fig. 3. Model SED corrected for interstellar extinction of 
A v = 12.8 and 2.13/im visibility function of WR 118 for 
the models Ml (solid lines), M2 (short-dashed lines), and M3 
(long-dashed lines) as given in Tabs. 0-g The squares refer to 
the observations (see text). 



and 14, respectively. The latter includes the data of Massey 
& Conti (1983) (V band), Cohen & Vogel (1978) (V band), 
Allen et al. (1977) (JHKL bands), Williams et al. (1987) 
JHKDM [8.4], [11.6], [12.5], [12.9], and [19] bands) and 
color-corrected IRAS photometry (IRAS Science Team 1988) 
([12] and [25] bands). Figures |5]andQ also show the 2.13 /im 
visibility and SED of those 5 models, which proved to be best 
suited for the dust shell of WR 1 1 8 in the framework of spher- 
ical symmetry. Their parameters are given in Tab. [j] and |[ For 
comparison, Williams et al. (1987) have found that the dust 
shell of WR 118 have the following parameters: r- m — 260i?*, 
y out = 30, T- m = 1410 K, p in = 4.9 • 10-19g/cm-3 and 
M dust =8.2- 10- 7 M Q /yr. 

The models Ml, M2, and M3 refer to an interstellar ex- 
tinction of Ay = 12.8, and the models Ml.l and M2.1 to 
Ay = 14.0. M3 is the model with a power law grain-size dis- 
tribution with 7 = 3. The other models refer to double-sized 
grains. In the latter case the relative abundance of 0.05 /im 
grains in number (no. 05) and mass (too. 05) are given in Tab. [j], 
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Fig. 4. Model SED corrected for interstellar extinction of 
Ay = 14.0 and 2.13 /im visibility function of WR 118 for 
the models Ml (solid lines), M2 (short-dashed lines), and M3 
(long-dashed lines) as given in Tabs. [||§ The squares refer to 
the observations (see text). 



which also contains the optical depth of the dust envelope Ty 
and its albedo at 2.13 /im (A^). 

Regarding the models with the double-sized grains (a = 
0.05 /im and a — 0.38 /im), the relative concentration of the 
0.05 /im grains is subject to some uncertainty, since the small 
grains determine the reddening, which is not precisely known. 
A more accurate determination of this parameter demands a 
more precise estimation of the optical brightness of WR 118 
at several wavelengths. However, the grain density at the inner 
boundary of the dust shell and the dust formation rate do not 
depend strongly on the adopted Ay and relative grain concen- 
tration. Thus, the estimations given in Tab. ^ can be considered 
as reliable in the framework of a spherically symmetric dust 
shell model. 

Furthermore, it turned out that a match of the SED at the 
region of 25 /im can be improved if the exponent in the power 
law of the dust density distribution is a = 2.1 instead of 2. It 
remains open if this deviation from the uniform outflow is sig- 
nificant with respect to the observational and theoretical uncer- 
tainties. The same effect as increasing a can be reached by de- 



creasing Y out to 300. This is the lower limit for the dust thick- 
ness above which the 25 /im flux is not affected any longer. 
However, such a sharp truncation of the outer parts of the dust 
shell looks less natural than, e.g., a gradual change of the mass- 
loss rate. 

The Ml and M3 models have almost the same albedo Afc, 
and the M2 and M3 models have almost the same circumstel- 
lar (est) reddening (Ay — Aj)/Ay\ cst . It defines the similarity 
of the Ml and M3 visibility curves and the similarity of corre- 
sponding SEDs (Fig. 1). The mass fraction of the grains with 
a > 0.3/rni in model M3 is ~ 50%. 

Table 3. Model dust-shell parameters for grain-size distribu- 
tions with n(a) oc a 1 and 0.005 /im< a < 0.6 /im. Ty is 
the optical depth at 0.55 /im; r- m , T ln and 9i n are the radius, 
the temperature (refering to the smallest grains) and the angu- 
lar diameter of the inner dust-shell rim, and Ak is the 2.2 /im 
dust-shell albedo. 
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Increasing the exponent 7 in the power-law grain-size dis- 
tribution from 3.0 (M3) to 3.5 (M4) leads to a decrease of 
the albedo A K from 0.36 (Tab. |) to 0.29, and to a signifi- 
cant increase of the visibility curvature, and thus to a worse 
match of the observations. The strong dependence of the visi- 
bility curve's shape on the albedo is confirmed by model cal- 
culations with double-sized grains. On the other hand, the de- 
crease of 7 to 2.5 does not change the albedo Ak much (to only 
~ 0.38) and correspondingly the visibility curve either (model 
M5). However, then the model V — J color decreases to 1.20 
(1.67 for M3) in contradiction to the observations. Moreover, 
the increase of the relative concentration of large grains en- 
forces a decrease of r m to keep a good match of the SED in the 
near-infrared. This, in turn, leads to a corresponding increase 
of temperature at the inner rim, T in , which is already very high 
(Tab.[|). Thus, from the models with power-law grain-size dis- 
tributions the best one is that with 7 = 3. 

If the wind of WR 1 18 is flattened due to some reason (due 
to stellar rotation, binarity), one can conclude from the circu- 
lar symmetric shape of the two-dimensional visibility function, 
that one looks face-on at the flattened shell ("disk"). Williams 
et al. (1987) estimated the fraction of WR 1 18's luminosity ab- 
sorbed by dust and re-radiated in the IR to be ~ 71%. If the 
star is not obscured by its dust disk, we estimate that this frac- 
tion should not be less than 50%. This means that the dust disk 
should be optically thick and cover more than one half of the 
sphere. Such a dust shell would not look like a geometrically 
thin disk, but rather as a sphere with a hole in the light of sight, 
i.e. as a torus with a large opening angle. However, first we 
should have clear indications of such asphericities, before we 
introduce a more complex geometry to the dust shell model of 
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WR 1 18. We note that van der Hucht (2001) suggested that all 
apparently single and persistently dust-making WC9 stars may 
possibly owe their heated circumstellar dust signatures to col- 
liding WC+OB wind effects. 

Recently, Veen et al. (1998) have attributed eclipse-like 
variations of the optical brightness of the dusty Wolf-Rayet 
stars WR103, WR121 and WR113 to the obscuration of 
starlight by dust clouds in the light-of-sight very close to the 
star, i.e. between the stellar surface and the permanent dust 
shell that is inferred from the IR excess. Moreover, they have 
estimated the sizes of the particles in the clouds to be of the 
order 0.1 ^tm. Crowther (1997) has discovered the visual fad- 
ing of WR 104, accompanied by the disappearance of high- 
ionization spectral lines, and placed the dust clouds inside the 
permanent dust shell. Thus, the decrease of the dust shell's in- 
ner boundary proposed by our modelling does not look as ex- 
traordinary as it appears with respect to the even more hostile 
conditions for grain condensation. 

5. Conclusions 

We presented the first diffraction-limited 2.13 [im speckle 
interferometric observations of the dusty Wolf-Rayet star 
WR 118 with 73 mas resolution. No evidence was found for 
any significant deviation of the dust shell from circular sym- 
metry. Radiative transfer calculations were carried out to model 
the spectral energy distribution and our 2.13 /im visibility func- 
tion. The main conclusion of these calculations is that the 
grains in the permanent dust shell of the persistently dust- 
making Wolf-Rayet star WR 118 grow to the relatively large 
sizes (a > 0.3 ^m). They contribute more than 50% in to- 
tal mass of the dust shell and rises its 2.13 ^m albedo up to 
~ 0.35. The observations can be fitted either with a model 
containing double-sized grains with radii of a = 0.05 /jm and 
0.38^m or, alternatively, by a model with a grain size distribu- 
tion function n(a) <~ a~ 3 , with a ranging between 0.005/im 
and 0.6^m. At the inner boundary of the dust shell, which 
has an angular diameter of 9i n = (17 ± l)mas and is lo- 
cated at r ln = (150 ± 30)i?„, the temperature of the small- 
est grains amounts to 1750 K ± 100K and the dust-shell den- 
sity to (1 ± 0.2) • 10~ 19 g/cm 3 . Adopting a wind velocity of 
Kvind = 1200 km/s, the dust formation rate is found to be 
(1.3 ±0.5) • 10- 7 M Q /yr. 

The existence of large grains in the dust shell leads to a 
revision of the previously assumed size of the inner dust-shell 
boundary. The inner dust-shell boundary is then located even 
closer to the star leading to higher grain temperatures. This, 
in turn, further worsens the conditions of grain condensation 
which are already extremely harsh in any case. The question 
how dust grains permanently condense under such conditions is 
still yet open (see Cherchneff et al. 2000 and references therein) 
but the existence of dust shells around Wolf-Rayet stars as in 
the case of WR 118 clearly shows its warrant urgency. 
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